General Relativistic Description of the Observed Galaxy Power Spectrum: 
Do We Understand What We Measure? 
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We extend the general relativistic description of galaxy clustering developed in Yoo, Fitzpatrick, and Zal- 
darriaga (2009). For the first time we provide a fully general relativistic description of the observed matter 
power spectrum and the observed galaxy power spectrum with the linear bias ansatz. It is significantly different 
from the standard Newtonian description on large scales and especially its measurements on large scales can be 
misinterpreted as the detection of the primordial non-Gaussianity even in the absence thereof. The key differ- 
ence in the observed galaxy power spectrum arises from the real-space matter fluctuation defined as the matter 
fluctuation at the hypersurface of the observed redshift. As opposed to the standard description, the shape of 
the observed galaxy power spectrum evolves in redshift, providing additional cosmological information. While 
the systematic errors in the standard Newtonian description are negligible in the current galaxy surveys at low 
redshift, correct general relativistic description is essential for understanding the galaxy power spectrum mea- 
surements on large scales in future surveys with redshift depth z > 3. We discuss ways to improve the detection 
significance in the current galaxy surveys and comment on applications of our general relativistic formalism in 
future surveys. 

PACS numbers: 98.80.-k,98.65.-r,98.80.Jk,98.62.Py 
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I. INTRODUCTION 

Measurements of the galaxy power spectrum can be used 
to infer the shape of the primordial matter power spectrum, 
providing valuable clues to understand the initial conditions 
of early universe. In the linear regime, galaxy bias — the re- 
lation between the galaxy and the underlying matter distribu- 
tions — is fairly generic and simple that the shape of the ob- 
served galaxy power spectrum reflects the shape of the matter 
power spectrum [1]. Therefore, recent theoretical work [2-5] 
has focused on interpreting the galaxy power spectrum mea- 
surements from quasilinear scales to nonlinear scales where 
the measurement precision is highest, yet galaxy bias is sig- 
nificantly affected by the complex nature of galaxy formation 
physics. 

However, the past decades have seen a rapid growth in this 
field, and large-scale galaxy surveys such as the Sloan Digi- 
tal Sky Survey (SDSS; |6[]) now cover a substantial fraction 
of the entire sky over a range of redshift, allowing for mea- 
surements of the galaxy power spectrum on sufficiently large 
scales with unprecedented precision J7H9|]. Recently, Dalai 
et al. II Ofl showed that the shape of the galaxy power spec- 
trum in the linear regime contains a very distinctive feature 
in the presence of the primordial non-Gaussianity. These 
very large-scale modes of the galaxy power spectrum are 
already measured with high signal-to-noise ratios from the 
SDSS, putting tight constraints on the amplitude of the pri- 
mordial non-Gaussianity [11], comparable to the limits from 
the Wilkinson Microwave Anisotropy Probe (WMAP). Fur- 
ther theoretical work ITll - [T7ll shows that future galaxy sur- 
veys can measure the galaxy power spectrum on large scales 
with the accuracy enough to detect the small non-Gaussianity 
from the simplest single field inflationary models [18] and the 
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nonlinear evolution of the Gaussian spectrum [ 19]. 

The contribution of the primordial non-Gaussianity to the 
galaxy power spectrum mainly arises from the gravitational 
potential and this relativistic contribution becomes compara- 
ble to the contribution of the matter fluctuation on large scales. 
However, at this large scale, where the relativistic effects be- 
come dominant, the standard Newtonian description of the 
galaxy power spectrum breaks down and the general relativis- 
tic description is therefore essential for understanding the ob- 
served galaxy power spectrum and deriving correct constraints 
from the measurements. Furthermore, there exists subtlety 
in theoretical calculations in this regime, where general rel- 
ativistic effects are important: Since the general covariance 
provides a large number of degrees of freedom, many theo- 
retical descriptions of observable quantities turn out to have 
unphysical gauge freedoms and they fail to correctly describe 
observable quantities that we measure on large scales. Natu- 
rally, these systematic errors in those theoretical calculations 
are negligible in the Newtonian limit, but they become sub- 
stantial on large scales 12011 . 

Fully general relativistic descriptions should be constructed 
by using observable quantities, rather than theoretically 
convenient but unobservable quantities as they are gauge- 
dependent. The gauge-invariance is a necessary condition for 
observable quantities and it should be used for the consistency 
check of theoretical calculations of observable quantities. For 
example, the observable quantities in the galaxy clustering 
case are the observed redshift, the observed galaxy position 
on the sky, and the total number of observed galaxies. Yoo, 
Fitzpatrick, and Zaldarriaga [20] developed a general rela- 
tivistic description of galaxy clustering and applied it to the 
cross-correlation of large-scale structure with CMB tempera- 
ture anisotropies, in which the largest scale modes can be most 
effectively probed. On low angular multipoles, the correct rel- 
ativistic prediction of the cross-correlation is larger by about 
a factor two than the standard Newtonian prediction, alleviat- 
ing the discrepancy between the theoretical prediction and the 
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anomalously large signal measured from the SDSS and the 
WMAP data (e.g., see HI). 

Here we compute the observed galaxy power spectrum, ac- 
counting for all the general relativistic effects but we only fo- 
cus on the linear regime, where the general relativistic effect is 
substantial. The observed galaxy power spectrum shows sig- 
nificant difference on large scales, compared to the standard 
description, i.e., the biased matter power spectrum in the syn- 
chronous gauge with the redshift-space distortion effect 12211 . 
The correct "real-space" matter power spectrum is anisotropic 
and the shape of the observed galaxy power spectrum changes 
with redshift. Especially, its measurements on large scales can 
result in a false detection of the primordial non-Gaussianity 
even in the absence thereof. 

The organization of this paper is as follows. In Sec. iHlwe 
compute the real-space matter power spectrum accounting for 
the relativistic effects and we discuss its stark difference from 
the matter power spectrum in the synchronous gauge. We 
extend the calculation of the real-space matter power spec- 
trum to computing the full observed galaxy power spectrum 
in Sec. MI Al and we forecast the detectability of the departure 
of the observed galaxy power spectrum from the standard pre- 
diction in the current and future galaxy surveys in Sec. lIIIBl 
We summarize our new findings and discuss the implication 
of our results for the proposed future surveys in Sec. lIVI 

The detailed and technical calculations of our results are 
presented in two Appendices. In Appendix lAl we summarize 
our notation convention and discuss the gauge issues associ- 
ated with metric representations. Calculations of observable 
quantities are presented in AppendixlBlwith particular empha- 
sis on the gauge-invariance of each equation. In computing 
our results, we assume there is no vector or tensor mode, while 
we present general formulas with vector and tensor modes. 
We use the Boltzmann code CMBFast 112311 to obtain the trans- 
fer functions of perturbation variables. 



II. REAL-SPACE MATTER POWER SPECTRUM 



bias factor b separating it from the rest of the contributions in 
Eq. (Q} and its simple but gauge-invariant structure provides a 
guideline for computing the full observed galaxy power spec- 
trum using Eq. (fTJ. More importantly, P ms ,(k) is the domi- 
nant contribution to the observed galaxy power spectrum on 
all scales. 

Therefore, we consider the matter fluctuation at the hyper- 
surface defined by the observed redshift z 

m dm n dm (z,n) = p m {t)(l + S m ) = p m (z)(l + m Sz ) , (2) 

where n is the observed position on the sky, Sz is the lapse 
in the observed redshift (see Appendix [Bj and we assume 
TOdm = 1 hereafter. Of particular importance is the distinction 
between S m and mg z - The former S m , commonly referred to 
as the matter fluctuation, is a gauge-dependent quantity and its 
meaning is ambiguous because the time slicing is unspecified. 
For example, the same matter density p m can result in differ- 
ent values of S m due to the change in the mean matter density 
p m (t) in different coordinate systems (note that p m (t) is just 
a function of coordinate time t). However, the latter ms z , 
also we call the matter fluctuation (at the observed redshift), 
is gauge-invariant and unambiguously defined, since the time 
slicing is set by the observed redshift z, rather than by an arbi- 
trary choice of coordinate systems as for S m . This perspective 
provides a theoretical ground for what ought to be the correct 
quantity for the so-called "real-space" matter power spectrum 

In computing the matter power spectrum, we need to com- 
pute the observed mean (ndm)n an d the observed matter fluc- 
tuation rhsz = n<im/(ftdm)n — I- 2 The nee d to distinguish 
rhgz from mg z arises from the ambiguity of quantities at ori- 
gin in msz- The full expression of m$ z is (see AppendixIBl 



(5 m + 3Hx)-3{-H6t + H x )o - 3 



3 / dr 



,(3) 



The general relativistic description of the observed galaxy 
number density n g is derived in 112011 and we rearrange the ex- 
pression in terms of gauge-invariant variables defined in Ap- 
pendixlAlas 



n p (z) 

+3 Sz x + 2 



1 + b m Sz + a x + 2 ip x + V - C afi e a e p (1) 
51Z „ d ( 5z y 



H— ( "~' x 
dz\H 



hp 6V L -2K, 



where n p is the mean galaxy number density in a homoge- 
neous universe and each term in the square bracket represents 
contributions from the volume distortion and the source effect 
(see AppendixIBl. 

In this section we compute the contribution of the matter 
fluctuation ms z = 8 m — 3 6z to the observed galaxy power 
spectrum. While the matter fluctuation itself is not directly 
observable, computing the matter power spectrum P mSz (k) 
merits close investigation, since ms z is boosted by a galaxy 



1 There are two widely known gauge-invariant variables for the matter fluc- 
tuation e m and e g defined in Bardeen 1 24] , where e m = S m + 3(1 + 
w)U{v + k/3)/k and e g = S m + 3(1 + w)H(P + 7') (see AppendixfA] 
for our notation convention). They describe the matter fluctuation in the 
matter rest frame (v + k/3 = 0) and in the zero shear frame (x = 0), 
respectively (e m and e 9 correspond to S v and 8 X in our notation). How- 
ever, in practice there is no theoretical reason to prefer one gauge condition 
to other conditions. The time slicing in observation is set by the (gauge- 
invariant) observed redshift z, and hence the matter fluctuation ms z in 
this frame correctly represents the "observed" matter fluctuation (note that 
Ttlgz = S m in the uniform-redshift gauge Sz = 0). Furthermore, all the 
above gauge-invariant variables e m , e g , and mg z are different on large 
scales (see Fig. [T}, and the gauge-invariance is not a sufficient condition 
for observable quantities. 

2 For an observed number density field n(z, n) at the observed redshift 2, the 
observed mean (n) n is obtained by averaging n over all angle n within the 
survey area Q, i.e., averaging over the hypersurface of simultaneity set by 
the observed redshift z. However, the usual ensemble average, also com- 
monly denoted as (n), is obtained by averaging n over the hypersurface set 
by the coordinate time t, and hence it is gauge-dependent. The observed 
mean and the theoretical mean (ensemble average) should be carefully dis- 
tinguished. 
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and it contains quantities evaluated at origin such as the time 
lapse St and so on. However, as noted in Yoo, Fitzpatrick, 
and Zaldarriaga [20], these quantities are nuisance in observa- 
tion due to their independence of the observed angle n: They 
are absorbed to the observed mean (ndm)n and subtracted 
from iidm to give rhs z without the ambiguities at origin (note 
that in observation rh$ z is derived from the observed number 
density nd m )- 

Though implicit in most theoretical calculations, we explic- 
itly account for this observational procedure. The observed 
mean number density is therefore 



/ d 2 nn d mO,n) / / d 2 n 
Jn I Jn 



(4) 



pm(z) 



1 - 3 {Ho 6t + H oXo ) + 3 (V - a x ) 



where V — V a e a . In the last equality, we assumed that per- 
turbation variables that vary with the observed angle n would 
vanish when averaged over angle within the survey area fi. 
However, this assumption is valid only if the survey area is 
infinite: Modes of wavelength larger than the survey scale 
would not average out but set a constant floor in the observed 
mean. Collectively, Eq. may have additional constant. 
However, since this additional constant is also independent of 
observed position n and gauge-invariant, it has no impact on 
the power spectrum computation, and we can safely assume 
that there is no further contribution in Eq. @. Also, note the 
subtle difference (ra dm ) n ^ Pm{z) ^ (n dm ) = p m (t), but 
their difference is at the level of metric perturbations. 

Therefore, the observed matter fluctuation can be written as 
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FIG. 1: Matter power spectrum Pm Sz (k, fit) at various redshifts. 
Upper panel: Solid lines represent the matter power spectrum at 
z = 0, computed by using Eq. Q along the line-of-sight direc- 
tion (pk = 1; thick) and along the transverse direction (fit ~ 0; 
thin). For reference, various lines indicated in the legend show power 
spectra of perturbation variables in the conformal Newtonian gauge 
and the synchronous gauge. Bottom panel: matter power spec- 
trum Pfi lgx (k, fik) at z > 0, but with its amplitude normalized to 
match Prh S z (k, /x^) at z = 0. Solid and dashed lines represent 
P-m 5z {k, Hk) with /ifc = 1 and /j,k = 0, respectively. The horizon 
scale at z — is shown as a vertical line. 



m Sz = n dm / (ridm)fi - 1 = (S m + 3#x) + 3 (a x - V) 



+ 3 / dr 
Jo 



K-^)'-(^ + Q^ 



(5) 



Now we can make a gauge choice and compute the matter 
power spectrum P 7 % gz (k). 3 By inspecting Eq. (0, it proves 
most convenient to compute the matter power spectrum in the 
conformal Newtonian gauge (\ — 0), 



8fL+3ib-3V + 3 



dr 



H afi e e H 



(6) 

where a x = a (= ip), cp x = ip (= (f>) in the conformal 
Newtonian gauge. We emphasize that mg z computed by us- 
ing Eq. (|6]l is gauge-invariant, i.e., its value is identical to those 
computed by using Eq. <(5j with any other choice of gauge con- 
ditions. Finally, the matter power spectrum can be computed 



3 Choosing a gauge condition per se has nothing unphysical nor incorrect. In 
fact, it is desirable for numerical convenience. The key point is that one has 
to be careful in relating theoretical predictions to observable quantities, for 
which the gauge-invariance is a necessary condition. However, by choos- 
ing a gauge condition before deriving equations, one loses an explicit way 
to check the gauge-invariance of theoretical predictions, and the predictions 
become gauge-dependent when they involve unobservable quantities. 



as 



■Ems* (k 



p (*)T*,.CM)i£ to (M) (7) 

~ P£(k) + 9 P^k) + 6 P m4 ,(k) + 9 n\ P v (k) , 

where P^(fc) is the power spectrum of the primordial curva- 
ture perturbations and Tm Sz is the transfer function for rhs z 
(we will use Tx to refer to the transfer function for the per- 
turbation variable X). We ignored the complication due to 
survey geometry, and the redshift dependence of each power 
spectrum term in the second line is suppressed for notational 
simplicity. Since projected quantities such as the integrated 
Sachs-Wolfe effect in Eq. ((6]) are intrinsically angular quanti- 
ties, they are nearly independent of the line-of-sight mode k" 
due to the Limber cancellation of the small scale modes, and 
we have ignored these contributions in computing Eq. (0. 

Solid lines in the upper panel of Fig. [T] represent the mat- 
ter power spectrum Pm Sz (k, A*fe) along the line-of-sight and 
the transverse directions. With the time slicing set by the 
observed redshift z, the "real-space" matter power spectrum 
P,n Sz (k) is no longer isotropic, and it is neither the matter 
power spectrum Pg(k) in the synchronous gauge (dotted), 
nor Pg (k) in the conformal Newtonian gauge (dot-dashed). 
Since the synchronous gauge corresponds to the dark matter 
rest frame, Pg (k) would be what we measure if we knew their 
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positions without measuring redshift z. However, the (coor- 
dinate) positions in the synchronous gauge are also gauge- 
dependent quantities, and we need the observed redshift z to 
define the time slicing and their observed positions. On small 
scales k ^> 10H, the matter power spectrum Pm Sz (k, fik) be- 
comes virtually isotropic, and the difference in Ps(k) with 
various gauge conditions vanishes compared to the correct 
matter power spectrum Pm Sz (k, 

The bottom panel shows the matter power spectrum 
Pm s ,(k, Hk) at high redshifts z = 0.5, 1, 2, 4.5, 8 along 
the line-of-sight direction (solid) and along the transverse di- 
rection (dashed). Compared to P/h Slt (k ) at z = 0, sig- 
nificant difference is immediately noticeable along the trans- 
verse direction, where the matter power spectrum has the zero- 
crossing scale at each redshift. Along the transverse direction, 
the matter fluctuation is mg z ~ + 3 -0. Since the mat- 
ter fluctuation 6% decreases with redshift but the gravitational 
potential if) (< 0) remains nearly constant, the matter fluctu- 
ation rhs z becomes zero and changes its sign on large scales, 
and this zero-crossing scale becomes smaller at higher red- 
shift. Physically, the magnitude of the lapse Sz in the ob- 
served redshift is larger than the density fluctuation 6 m on 
large scales, such that the change in the mean matter density 
pm (t) overwhelms the density fluctuation and hence m$ z ap- 
pears to be underdense in observation. The matter power spec- 
trum Pm 5z {k" ) along the line-of-sight direction maintains the 
similar trend at z = shown in the upper panel. However, 
compared to the matter power spectrum Pg(k) in the syn- 
chronous gauge (dotted), the turn-off scale of Pfa 5z (k,Hk) be- 
comes smaller, as the horizon size decreases with redshift. 



III. OBSERVED GALAXY POWER SPECTRUM 

A. Theoretical Prediction 

Having computed the real-space matter power spectrum in 
Sec. lU we are now in a position to tackle the full complexity 
of the observed galaxy number density n g in Eq. ((T). Fol- 
lowing the same procedure, we compute the observed mean 
galaxy number density by considering contributions indepen- 
dent of the observed galaxy position, 



ln g )n = n. 



l+b ms zo +'d Sz Xo +2 -—hp 5V Lo -2 JC a 



where the surface terms in each component are 

5z Xo = (H 5t + HoXo) - (V - a x ) Q , 
mszo = -3 Sz Xo , 

STZ = (xo + St ) - 6Zx ° 



(8) 
(9) 



H ' 

SV Lo = (H St + H oX o) + 



Xo + 5t 



Therefore, the observed galaxy fluctuation can be obtained by 
subtracting the observed mean as 

&g = n g /(n g ) n - 1 = b (m Sz ~ m Szo ) + a x + 2 p x (10) 



+V - C af} e a e p + 3 (5z x - 5z Xo ) + 2 



-H-j^ - hp (SV L - 5V Lo ) -2{K-K ). 

This equation maintains the gauge-invariant structure, yet 
lacks any ambiguity at origin. 

As discussed in Sec. [II] we make simplifications to facilitate 
computing the observed galaxy power spectrum P g (k, (ik) by 
adopting a simple survey geometry, ignoring contributions 
from the projected quantities, and choosing the conformal 
Newtonian gauge (x = 0). The observed galaxy power spec- 
trum is therefore 



P g (k)=P (fc)T g (k,z)T g *(k,z), 
with the transfer function 



T„ 



b (T„; + 3 
+3 (T y - T» 



3 TV) + T* + 2 T t 



-hp \ T V - P4, 
bT" +NT^ 



Ttjj — Ty 

Ur s 

Tip — Ty 



-4 

dz 



(11) 



Ty — T^ 
H 



kT, 



'-■4, 



ip, k MT v 



where two coefficients in the transfer function are 

1 + z dH . ,„ „ s (hp 



(12) 



7^ = 36-1 



M = 3b 



1 



H dz 
zdH 



(hp - 2) - 



(hp - 2) - 



Hr. 
hp -2 



H dz ' ' V Ur. 

Finally, the observed galaxy power spectrum in Eq. ( fTTT i is 

P g (k,fi k ) = b 2 P^(k)+N 2 P i ,(k) + 2bNP m i P (k) 



(13) 



/4 



M 2 P v (k) ^P mv (k) - ^P^(fc) 



n 



n 



4 ^Pv(k) , 



(14) 



where we ignored the contribution from Tp /U, which van- 
ishes quickly at z > 0. With our fiducial cosmology these co- 
efficients range M - 3 b = -0.5 - 0.5 and M - 3 b = 0.8 ~ 
1.5 for the luminosity function slope hp — 2 of the source 
galaxy sample, where the source effect is minimized. 4 The ob- 
served galaxy power spectrum P g (k, fik) is therefore approxi- 
mately the biased matter power spectrum b 2 Pm Sz (k, Hk) with 



4 The source effect cannot be exactly cancelled on all scales, even for the 
case 5p = 2, because the fluctuation 8T>l in the luminosity distance 
reduces to the standard convergence k only in the Newtonian limit, i.e., 
5p ST>l + 2 K 7^ on large scales, even when bp = 2. 
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FIG. 2: Observed galaxy power spectrum P g (k, fj,k) atz = 1. Thick 
lines represent the observed galaxy power spectrum, computed by 
using the full general relativistic description in Eq. dl4b along the 
line-of-sight direction (pk = 1; solid) and along the transverse di- 
rection (/ife = 0; dashed). Dotted lines show the observed galaxy 
power spectrum using the standard method in Eq. d 1 5 b . For compari- 
son, the matter power spectrum 6 2 Pm 5z (k, /j,k) is shown as thin lines 
with the galaxy bias factor 6 = 2. 



the standard redshift-space distortion factor, in which case two 
coefficients are exactly J\f = 3 b, Ai ~ 3 b (when 5p = 2 is 
assumed). 

In comparison, the observed galaxy power spectrum and its 
transfer function in the standard method are 




aj - 



0.001 0.01 
wavenumber {hMpc" 1 ) 

FIG. 3: Systematic errors in theoretical modeling of the observed 
galaxy power spectrum. Upper panel: fractional difference of 
P g (k, fj,k) at various redshift slices, compared to the standard de- 
scription P st d{k, Hk) in Eq. M5\ along the line-of-sight direction 
(solid) and along the transverse direction (dashed). Bottom panel: 
detection significance of the departure from the standard descrip- 
tion in a cosmic-variance limited survey as a function of maximum 
wavenumber. The survey volume is divided into four spherical shells 
with redshift range z = ~ 1 (bottom solid), 1 ~ 3, 3 ~ 6, and 
6 ~ 10 (top solid), and the signal is computed at the mean redshift 
of each shell. 



P std (Mfc) = P*{k)T std (k,z)T* t<i (k,z) (15) 
= b 2 P s s (k) - A 2 -^P Sv {k) + 4 ^P v (k) , 

and 

kT 

T atd = bT%- f 4^ , (16) 

respectively. The standard transfer function T stc \ is composed 
of two transfer functions (T r f t and T v ) in two different gauge 
conditions, one in the synchronous gauge and one in the con- 
formal Newtonian gauge. 

Figure [2] shows the observed galaxy power spectrum 
P g (k,Hk) at z = 1, computed by using Eq. ( fT4l i. As we 
noted in Sec. [TTJ the matter power spectrum Pm Sz (k, /ife) is 
the dominant contribution to the observed galaxy power spec- 
trum on all scales, and the redshift-space distortion from 
the volume distortion enhances the power. Along the trans- 
verse direction (thick dashed), the observed galaxy power 
spectrum P g (k ± ) largely traces the biased real-space matter 
power spectrum b 2 P 7 ~ nSr (k ) (thin dashed) as the standard de- 
scription of the observed galaxy power spectrum (lower dot- 
ted) is the biased matter power spectrum in the synchronous 
gauge P s td{k ± ) = b 2 Pf(k). However, on large scales, 
the standard description overestimates the observed galaxy 
power spectrum along the transverse direction due to the 
contribution of 8z in Pm Sll {k ± )- Similarly, the observed 



galaxy power spectrum P g (k") along the line-of-sight direc- 
tion (thick solid) significantly deviates from the standard de- 
scription P s td(fc") = (b + f) 2 Ps(k) (upper dotted) on large 
scales, where / is the growth rate of structure. The deviation 
arises mainly from the difference between b 2 P 7 - nSr (k") and 
b 2 Pg(k), but there exists small contribution from the volume 
distortion in addition to the standard redshift-space distortion 
effect. 



Figure [3^ compares the correct galaxy power spectrum 
Pg(k, Hk) in Eq. (fT4b with respect to the standard description 
Pstd(k,fJ,k) in Eq. ( fTBT ) at various redshift slices. Solid and 
dashed lines represent the fractional difference AP g /P st d — 
Pg/Pstd — 1 along the line-of-sight direction (/ife = 1; solid) 
and along the transverse direction (/ife = 0; dashed), respec- 
tively. On large scales k < 0.01/iMpc -1 , the standard de- 
scription underestimates the observed galaxy power spectrum 
along the line-of-sight direction and overestimates along the 
transverse direction. At higher redshift, the fractional dif- 
ference in the observed galaxy power spectrum at a fixed 
wavenumber is progressively amplified, because the relativis- 
tic effects increase as the horizon size becomes smaller. A 
few percent deviation is already present at k — O.Ol/iMpc -1 
at z > 4.5. 
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B. Statistical uncertainties on power spectrum measurements 

Of the utmost importance is the detectability of the sys- 
tematic errors in theoretical modeling of the observed galaxy 
power spectrum P g (k,Hk) on large scales. We consider sta- 
tistical uncertainties on measurements of the two-dimensional 
anisotropic power spectrum by using a simple mode counting 
approximation. Since Fourier modes are effectively indepen- 
dent of each other, the statistical uncertainty decreases with 
1 / y/Nk, where Nf. is the number of Fourier modes available 
given a survey volume V s . However, sampling of individual 
galaxies contributes the shot noise errors l/(n s )n to the sig- 
nal P g (k, /ife), and the number of independent Fourier modes 
is in fact a half of Nf. as 5 g (k) represents a real quantity in 
configuration space. Therefore, the statistical uncertainty of 
each Fourier mode is 11251 uol 



(17) 



V k V cS 



and the detection significance of the departure in the observed 
galaxy power spectrum from the standard method prediction 
is 



l\ 2 = (AP g " 2 
N 



where AP = P„ - P atd , 



27rA(fc- L ) 2 AfcH 
(2nf 



(19) 



is the Fourier volume accounting for the positive and negative 
k", and the effective survey volume is 



V cS 



dV s 



( n g)n P g 



1 + (n g )nP g 



(20) 



For definiteness we consider galaxy samples mea- 
sured from a large-scale survey with full sky coverage 
(cosmic-variance limited) as a function of survey redshift 
depth. We choose the galaxy number density (n 9 )n = 
10 -4 (/i -1 Mpc)~ 3 , similar to what the current large-scale sur- 
veys aim to measure (e.g., (2j\\ ). Since the shot noise contri- 
bution is small (n 9 )nP 9 > 1 on scales of our interest k < 
0.01/iMpc~ , the effective volume in Eq. (|20T > is therefore 
identical to the survey volume V e s — V s . Here we assume 
a fiat ACDM universe with the matter density Q, m = 0.24, the 
baryon density fib = 0.042, the Hubble constant h = 0.73, 
the spectral index n s = 0.954, the optical depth to the last 
scattering surface s — 0.09, and the primordial curvature per- 
turbation amplitude = 2.38 x 10~ 9 at k = 0.05 Mpc" 1 
(a 8 = 0.81). 

The significance (S/N) 2 of the systematic errors as a func- 
tion of maximum wavenumber is shown in Fig. [3J3, where we 
divide the survey volume into four spherical shells, covering 
the full sky with redshift range z = 0~l,l~3, 3 ~ 6, and 
6 ~ 10. At the lowest redshift bin (z — ~ 1), the survey 



volume is V s — 58 (h 1 Gpc) 3 and the minimum wavenum- 
ber is fc min = Ak = 2tt/V s 1/3 = 0.0016/iMpcT 1 , at which 
the fractional difference is 8% along the line-of-sight direc- 
tion (solid) and is close to zero along the transverse direction 
(dashed) as shown in Fig.[3k. With only a handful of modes at 
fc m in, the systematic errors in theoretical modeling of the ob- 
served galaxy power spectrum is negligible at z < 1 (bottom 
solid line in Fig. |3p). Similar conclusion can be derived for 
the survey volume at z = 1 ~ 3 (second solid line from the 
bottom). 

However, at the higher redshift bins (two upper solid lines), 
there exist two critical differences in assessing the significance 
of the systematic errors: Larger survey volume and smaller 
horizon size. About a factor of ten larger volume is avail- 
able in the higher redshift bins than in the lowest redshift 
bin, which reduces the statistical uncertainties at a fixed fc by 
10 1/2 ~ 3.2 and provides fc min smaller by 10 1/3 ~ 2.2. Fur- 
thermore, the relativistic effects are larger at a fixed k due to 
the smaller horizon size at the higher redshift bins, increasing 
the fractional difference AP g /P st d at each fc. Consequently, 
the significance of the systematic errors increases dramati- 
cally with survey redshift depth. 

In conclusion, the volume available at z < 3 is less than 
the Hubble volume and hence the standard Newtonian de- 
scription is statistically indistinguishable from the general rel- 
ativistic description of the observed galaxy power spectrum at 
z < 3. However, the systematic errors in theoretical modeling 
increase substantially with redshift and correct general rela- 
tivistic description is essential in understanding the observed 
galaxy power spectrum at z > 3. 

In practice, further complication will be present for mea- 
suring the largest scale modes in galaxy surveys, because they 
are affected by survey window functions and large contiguous 
region is required. Long line-of-sight modes are relatively im- 
mune to these difficulties, while modeling the redshift evolu- 
tion of galaxy bias may complicate the theoretical interpreta- 
tion of the observed galaxy power spectrum P g (feH) along the 
line-of-sight direction. Proper application of our method to 
observations will require further investigation with close ties 
to the specifications of a given survey geometry. 



IV. DISCUSSION 

We have extended the general relativistic description of 
galaxy clustering developed in Yoo, Fitzpatrick, and Zaldar- 
riaga B20I1 and we have provided, for the first time, the correct 
general relativistic description of the observed galaxy power 
spectrum with particular emphasis on the expressions for ob- 
servable quantities and their gauge-invariance. In the galaxy 
clustering case, the observable quantities involve the observed 
redshift z, the galaxy position n on the sky, the number 7Y tot 
of observed galaxies, and its apparent luminosity L. These 
quantities are different from the quantities defined in a ho- 
mogeneous and isotropic universe, which are theoretically 
more tractable, yet unobservable (gauge-dependent) quanti- 
ties. The gauge-invariance is a necessary condition for ob- 
servable quantities and the observed galaxy power spectrum 
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should be expressed in terms of the observable quantities. 

There are two key contributions to the observed galaxy 
power spectrum: The real-space matter fluctuation and the 
volume distortion. While the standard redshift-space distor- 
tion effect accounts for most of the contributions from the vol- 
ume distortion, the former requires a careful definition of its 
meaning and it exhibits significant difference on large scales 
compared to the standard real-space matter fluctuation. Since 
the observed redshift z defines the hypersurface of simultane- 
ity in observation, the correct description of the real-space 
matter fluctuation is the matter fluctuation at the observed 
redshift, mg z = S m — 3 8z (mg z — S m in the uniform- 
redshift frame 5z = 0). As it is defined by observable quanti- 
ties, msz is naturally gauge-invariant, in contrast to the usual 
(gauge-dependent) matter fluctuation S m , which requires a 
specification of gauge conditions or coordinate systems and 
differs in its value contingent upon the gauge choice. 

As the observed redshift z depends on angle, the real-space 
matter power spectrum P mSz (k) is no longer isotropic: Com- 
pared to the matter power spectrum Pg (k) in the synchronous 
gauge, P m5z (k) is enhanced along the line-of-sight direction 
but it is suppressed along the transverse direction, because the 
lapse Sz in the observed redshift can increase or decrease the 
mean matter density, affecting S m in both direction. Conse- 
quently, the shapes of the real-space matter power spectrum 
Pm 5z (k) and the observed galaxy power spectrum vary with 
redshift at k < 0.01/iMpc , while the shape of the matter 
power spectrum Pg(k) in the synchronous gauge is indepen- 
dent of redshift. Therefore, a proper modification is necessary 
for methods to use the multipole moments of the observed 
galaxy power spectrum on large scales to constrain the growth 
factor of the matter fluctuation. Traditionally, measurements 
of the observed galaxy power spectrum at different redshift 
slices are combined to reduce the statistical uncertainties, but 
with the redshift evolution of the power spectrum shape on 
large scales, they can be used to extract additional cosmolog- 
ical information. 

Recently, it is shown ifioll that the primordial non- 
Gaussianity in gravitational potential can give rise to the scale- 
dependence of the observed galaxy power spectrum on large 
scales, where galaxy bias is assumed to be linear. Further- 
more, it is guaranteed lfl2l IyA l28ll that future galaxy sur- 
veys will detect non-Gaussian signatures on large scales, aris- 
ing from either primordial or nonlinear evolution. However, 
its signature on large scales is similar to the rising power 
of the real-space matter power spectrum P mSz (k) in Fig. Q] 
and the observed galaxy power spectrum in Fig. [2] since they 
all originate from the contribution of the gravitational poten- 
tial. Therefore, without proper theoretical modeling of the 
observed galaxy power spectrum, its measurements on large 
scales can be misinterpreted as the detection of the primor- 
dial non-Gaussianity even in the absence thereof, and its false 
detection can be used as the evidence against the simplest sin- 
gle field inflationary models. With the full general relativistic 
description of the observed galaxy power spectrum developed 



here, the effect of the primordial non-Gaussianity on the ob- 
served galaxy power spectrum will be investigated in a sepa- 
rate work. 

While the current galaxy surveys cover a large fraction of 
the entire sky, it is still difficult to directly measure the depar- 
ture of the observed galaxy power spectrum from the standard 
description at z < 3 with high significance, simply due to the 
limited number of large-scale modes. However, recent work 
ll29ll suggested a method to eliminate the cosmic variance er- 
rors by comparing two different biased tracers. The scale- 
dependence of the ratio of two galaxy power spectra arises 
from the volume distortion: Unlike in the standard description 
of the primordial non-Gaussianity, all biased traces depend on 
the same matter fluctuation mj z and hence there would be 
no scale-dependence, if it were not for the volume distortion. 
With substantial increase in the signal-to-noise ratio, current 
galaxy surveys may be able to detect the general relativistic 
effects in the observed galaxy power spectrum at low redshift. 

In future galaxy surveys with larger redshift depth, cor- 
rect general relativistic description should be an essential el- 
ement in power spectrum analysis, as significant systematic 
errors will be present in the standard description of the ob- 
served galaxy power spectrum due to smaller horizon size 
and larger survey volume at high redshift. Furthermore, our 
formalism can be easily extended to any alternative theories 
of gravity. Since these theories are identical in the Newto- 
nian limit but have distinctive relativistic effects, alternative 
theories of gravity can be most effectively probed on cosmo- 
logical scales. Comparison of their theoretical predictions to 
the galaxy power spectrum measurements on large scales can, 
therefore, provide a way to test general relativity and put con- 
straints on the alternative. 

While we considered a spectroscopic galaxy survey here, 
a majority of planned large-scale surveys such as the Dark 
Energy Survey (DES), the Panoramic Survey Telescope & 
Rapid Response System (PanSTARRS), and the Large Syn- 
optic Survey Telescope (LSST) lack of observationally ex- 
pensive spectroscopic instruments, and instead they rely on 
photometric redshift measurements. However, the degrada- 
tion of redshift measurement precision will have little impact 
on the detectability of large-scale modes, since it only results 
in radial smearing of galaxy positions (Ar < lOO/i Mpc at 
z = 1, and smaller at higher redshift) and the general rela- 
tivistic effect appears on large scales k < 0.01/iMpc - . 
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Appendix A: FLRW metric and notation convention 

We summarize our notation convention for the most general 
metric in an inhomogeneous Friedmann-Lemaitre-Robertson- 
Walker (FLRW) universe and we discuss the gauge trans- 
formation properties and gauge freedoms of the perturbation 
variables used in this paper. 

The background universe is described by a spatially homo- 
geneous and isotropic metric with scale factor a(t), 

ds 2 = g ab dx a dx b = -dt 2 + a 2 (t) g aj} dx a dx? , (Al) 

where g a p is the metric tensor for a three-space with a con- 
stant spatial curvature K = —H 2 (1 — f^tot)- We extensively 
use the conformal time r, defined as a(r) dr — dt . Through- 
out the paper, we use the Latin indices for the spacetime com- 
ponent and the Greek indices for the spatial component, and 
we set the speed of light c = 1 . 

Small departure of the metric tensor in the inhomogeneous 
universe from the background metric can be represented as a 
set of metric perturbations: 

5g 0Q = -2 a 2 a, (A2) 
Sgoa = -a 2 B a = -a 2 (/3. a + B a ) , 

5g a /3 = 2 a 2 C a/3 = 2 a 2 (ip g a p + 7 ^ + C( a \/3) + C a p) , 

where the vertical bar represents the covariant derivative with 
respect to spatial metric g a p , and the perturbation classifi- 
cation from its transformation properties is represented with 
single or double indicies for vector-type or tensor-type, re- 
spectively. We adopt the classification scheme introduced in 

The general covariance in general relativity guarantees that 
any coordinate system can be used to describe physical sys- 
tems and the choice of a coordinate system bears no physi- 
cal significance. A coordinate transformation accompanies a 
change in the correspondence of the real physical (inhomo- 
geneous) universe to the homogeneous background universe, 
known as the gauge transformation. Here we consider the 
most general coordinate transformation 

£ a = x a + C , (A3) 

where f° = (T, C a ) and C a = L' a + L a . Under the coordi- 
nate transformation, the scalar metric perturbations transform 

as 

a = a-T'-HT, (A4) 
P = 0-T + L', 
<fi = (p — HT, 
7 = 7-i , 
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where the conformal Hubble parameter is H — a' /a = aH 
and we denote the conformal time derivative as a prime. The 
vector metric perturbations transform as 



B a — B a + L' 



(A5) 



Since tensor harmonics are independent of tensors that can be 
constructed from coordinate transformations, tensor-type per- 
turbations remain unchanged under the gauge transformation. 

Based on the gauge transformation properties, fully gauge- 
invariant quantities can be constructed as 



a x 

fx 
v 



V 



a - x'/a , 

f - Hx , 

v + k (3 — k x/a 



(A6) 



for scalar perturbations, where x = a (P + l') and x = X ~ 
aT . Our notation is chosen to facilitate the choice of gauge 
conditions [30], e.g., a x — a, cp x = tp, and v x = v in the 
conformal Newtonian gauge (fi = 7 = 0) or the zero-shear 
gauge (x = 0). For vector perturbations, two gauge-invariant 
quantities are 



* Q = B a + C' a , 



(A7) 



and we have used V a = v x , a + V a in the text. These gauge- 
invariant variables correspond to <&a, &h, vf \ ^P, and v c in 
Bardeen's notation 12411 . In addition to these gauge-invariant 
variables, one can construct other gauge-invariant variables, 
e.g., S v = 5 + 3(1 + w)T~L(v + kf3) /k, where 5 V = S in the co- 
moving gauge, and indeed there are as many possibilities for 
gauge-invariant variables as for the choice of gauge condition. 

Due to the spatial homogeneity of the background universe, 
all the perturbations should be invariant under pure spatial 
gauge transformations (T — 0, C a ^ 0). Therefore, those 
perturbation variables (f3,-f,B a ,C a ) that transform with L 
or L a carry unphysical gauge freedoms and they can ap- 
pear in physical quantities only through the combinations x 
and \& Q that are invariant under spatial gauge transformations 
JUl]. When the observed mean is subtracted from the ob- 
served number density field in Sees. [IT] and [III] these pertur- 
bation variables (f3, 7, B a , C a ) may leave unphysical gauge 
freedoms in the observed fluctuation field ruining its gauge- 
invariance. Fortunately, the damage is relatively innocuous, 
affecting only the monopole and the dipole at origin. In the 
paper, we explicitly resolve this issue by rearranging physi- 
cal quantities in terms of spatially gauge-invariant combina- 
tions x an d 'fn: before any operation is taken. 



Appendix B: Derivation of the Gauge-Invariant Equations 

Here we derive the gauge-invariant equations used in the 
power spectrum analysis. The derivation closely follows 
our previous work (20J, but we pay particular attention to 



the gauge-invariance of perturbation variables of the observ- 
able quantities. Unphysical gauge terms are isolated and re- 
moved by expressing observable quantities in terms of gauge- 
invariant variables defined in AppendixlAl 

We parametrize the photon geodesic x a (v) in terms of the 
affine parameter v, and its wavevector is 



dr a 
dv 



-(1 + H ~-(e a +Se a ) 
a a 



(Bl) 



where v and e a represent the photon frequency and the pho- 
ton propagation direction measured from the observer in the 
homogeneous universe. The null condition (k a k a = 0) for 
the wavevector constrains the zeroth order propagation direc- 
tion e a e a = 1, and the geodesic equation (k b k a b — 0) yields 

v oc 1/a and e' a = e^e"^ . 

The dimensionless perturbations to the wavevector k a (v) 
are represented as 5v and Se a for each component of Eq. ( IB II ). 
and for the coordinate transformation in Eq. (|A3l l these pertur- 
bation components of the wavevector transform as 



S V = SV+ — T + ZHT, 
dX 



(B2) 



Se a = Se a + 2UTe a - Te' a - — C a ~ e a g C [i , 

dX 

where d/dX = d T — e a d a and it is related to the zeroth or- 
der photon path d/dv = k a d a = (v/a)(d/dX) . Based on 
the gauge transformation properties, we define two gauge- 
invariant variables for the wavevector as 



5v x = 5v 



(B3) 



be. 



5e a + 2H X e a 



a / dX ' H 



where Q a = 7'" + C a is a pure gauge term, transforming as 
g a = g a - C a . 

The dimensionless perturbations 8v and Se a to the 
wavevector are also subject to the null condition and satisfy 
the geodesic equation. In terms of the gauge-invariant vari- 
ables, the null condition is 



x 



fx - ^ a e c 



C a0 e a e 



(B4) 



and the temporal and spatial components of the geodesic equa- 
tion are 



dX 
and 



(6v x + 2a x ) = (a x - <p x )' - (* a]B + C' aB ) e a e? 



(B5) 
(B6) 



(5e^ +2ip x e a +^ a + 2 C%e p ) 

-e< (5e% + 2 y x e a + + 2 C^e' 3 ) 

= Seteg - 5v x e' a + [a x - <^ x ) |q - ^p la e p - C^ePe 1 , 



respectively. Fictitious gauge freedoms in 8v and Se a are 
completely removed, and Eqs. ( IB4I ). dB5l ), and JB6I ) are mani- 
festly gauge-invariant. 
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As our observable quantities are obtained by measuring 
photons from galaxies, we are mainly interested in perturba- 
tions along the photon geodesic given observed redshift z and 
observed angle n. Therefore, when we consider a coordinate 
transformation with the observable quantities fixed, the affine 
parameter A is also affected by the coordinate transformation 
in Eq. JA3I ), i.e., 



x a {X) = z a (A)+£ a (A) , 
and the deviation of the photon geodesic is 



8x a = Sx a + C ~ k a - SX , 
v 



(B7) 



(B8) 



where SX = A — A and dX = dX (1 — 2HT). The time lapse 
and the spatial shift are 



St q = 5t - 
Sx» = 5x a Q 

at the origin A = A c and 



T , 



(B9) 



6t s = St s + T s - SX S , 
Sx§ = 5x a s + L<* + e a SX s , 



at the source A = A s , where the subscripts represent that the 
quantities are evaluated at the origin x a (A ) or the source po- 
sition x a (X s ). 

A comoving observer with the four velocity u a = [(1 — 
a) /a, v a /a] measures the redshift parameter z from galaxies 
and it is related to the photon wavevector as 



_ [k a (X s )u a (X s )} _ 1 
+ [k*(\o)u a (\ )} -a s {L + 0Z) 

' '1 + HoSto 



(Bll) 



V-H X 



— {1-Hx + {H oX o + Uo5t ) 



dr 



V-a y 



(a x - <p x )' - {* a{ p + C a g) e a e? 



where r s is the comoving line-of-sight distance to the source, 
V = V a e a is the line-of-sight velocity, and a s — a[r(A s )]. 
The redshift parameter z^ of the source, commonly defined 
as 1 + Zh = l/a s , needs to be distinguished from the ob- 
served redshift z, and Eq. ( IB lib defines the lapse term 5z in 
the observed redshift. Spurious spatial gauge freedom in 6z 
is removed and its temporal gauge dependence is isolated 
Sz = Sz + HT. We define a gauge-invariant variable for 
the lapse in the observed redshift as Sz x = 5z + H\ ■ 

The spatial displacement of the photon path can be obtained 
by integrating the spatial components of the geodesic equa- 
tion in Eq. (IB6I > and using Eqs. dB9l ) and ( IB 10b . The radial 
displacement is 



Sr 



c a x s 



(B12) 



(Xo + 5t ) - ^ 
(Xo + 5t ) - ^ 



e a (Sx«-G a ) 
e a G a 



dX 5i 



dr (a x -<p x - * Q e" - C a/3 e a e p ) 



where (d/dX) Sx" = — S e". W e have used the null condition 
for the wavevector in Eq. ( IB4I ) and defined the radial displace- 
ment Sr with respect to the comoving line-of-sight distance 
r 8 , slightly different from the definition used in 12011 . Simi- 
larly, the angular displacements are 



= v n x„ = 



U^ x -Q a )\ o 



(B13) 



9 a G a 



dr 



8 



88 



(BIO) and 



r s 8<j> = cj) a x" = 



dr L{^ a +2C$e f) ) 



(B14) 



r sin ( 



8 



(a x -<p x - * a e Q - C a pe a eP) 



with the orthonormal direction vectors 8 and cf) perpendicular 
to the observed direction n. The gravitational lensing conver- 
gence is then defined as 



cot0 



88 



8 



e Q (<5e" + * a + 2C^) - 



(B15) 



dr - 



2r, 



dr 



8 4>a 8 



89 sind. 



(<f a + 2C2e /3 ) 



dr 



e a G a 



2 rr, 



V 2 (a x -<p x -* a e a -C afi e a ef > ) 



8 



b a 8 



a 88 sind. 



where the Laplacian operator on a unit sphere is V 2 = 
(8 2 /88 2 ) + cot 0(0/90) + (1/ sin 2 8) {8 2 /8<j) 2 ). The pres- 
ence of the gauge terms Q a in Eqs. dB12b -( lBT31 > explicitly 
indicates that the spatial displacements Sr, 88, and 8(f) and 
the convergence n are gauge-dependent, which can be un- 
derstood as they represent the difference between the ob- 
served source position and the unobservable unlensed source 
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position and subsequently the function thereof. For nota- 
tional convenience, we define two gauge-invariant quantities: 
SK = Sr + [e a g a ] s and K = K + [e a G a } s /r s - [6 a (d/d6) + 
($ a / S in9)(d/d<f>)](g°'/2). 

Finally, the volume occupied by the observed galaxies can 
be obtained by tracing backward the photon geodesic x a (X) 
given the observed redshift z and angle n, and the total num- 
ber of galaxies in the observed volume is 



— a dx a dx b dx c , 

-g n p Eabcd u — — — - dz dO , 

oz oO dm 



{l + z) 3 H 



dz d9 dcj) x 



-A7 + v a e a + '3Sz + 2 



Sr 



l + 3(p 

H^-Sr 

oz 



(Bi6) 



where n p is the physical number density of the source galax- 
ies, the metric determinant is y/—g = a 4 (1 + a + 3 ip + A7), 
and £ a bcd is the Levi-Civita symbol. Equation ( IB16b means 
that iVtot is rip times the physical volume described by the 
observed redshift z and angle n = (9, <fi), and it defines the 
observed galaxy number density n g as 



{1 + zfH 



dz dO 1 



(B17) 



providing the relation to the physical number density n p . 

With the contributions from the volume distortion between 
the observed and the physical, the observed galaxy number 
density has additional contributions from the physical galaxy 
number density, since more galaxies are observed along the 
overdense region due to the magnification given the threshold 
Fthr for the observed flux, and this source effect depends on 
the intrinsic luminosity function of the source galaxy popula- 



tion, drip/dL oc L 



n p (ithr)(l - 5 p SV L ), where 



the inferred luminosity threshold for the observed galaxies is 
Lthr = 47rZ?!(z)F t hr, the luminosity function slope in mag- 
nitude is p = 0.4 (s — 1), and the luminosity distance in a 
homogeneous and isotropic universe is Dl. 

Given the observed redshift z and angle n, the fluctuation 
5T>l in the luminosity distance T>l, defined as T>l(z,ti) = 
D L (z)(l + SV L )AsMM 



SV L 



UoSt 

1 
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and it is related to the magnification fj, as p, = 1 — 2 SV^. The 
last term in Eq. ( IB 1 8b is the source term for the expansion of 
the wavevector d — k a a in the conformally transformed met- 
ric g a b = {p/a)g a b with the corresponding affine parameter A, 



8(R ab k a k b ) = A(a x - <p x ) - a xMp e a e 13 - 2 

+ (K\^ - A* Q )e" + ((%p - AC a p)e a e? , (B19) 



where R a b and A are the Ricci tensor and Laplacian operator 
of a three-space. 

The last ingredient of our formalism is the linear bias ap- 
proximation that relates the physical galaxy number density 
rip to the underlying matter distribution p m . Observation- 
ally, the time slicing is set by the observed redshift z and this 
choice is the only physical way in cosmological observations 
for defining the hypersurface of simultaneity without gauge 
ambiguity. The matter density at the observed redshift z of 
source galaxies is p m = p m (t)(l + 8 m ) = p m (z)(l + m Sz ) 
where m$ z = 8 m — 3 Sz and it is simply the matter fluc- 
tuation m$ z = S m in the uniform-redshift gauge (Sz = 0). 
Therefore, the simplest linear bias ansatz we adopt is n p = 
n p (z)(l + b m Sz ) . 

Finally the observed galaxy number density can be written 
in a manifestly gauge-invariant manner as 



lg = Up(z) 



1 + 6 msz + a x + 2 f x + V 
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and this completes the derivation of Eq. (d}. In Eq. ( |B20| > 
n p (z) is the mean galaxy number density in a homogeneous 
and isotropic universe, in a close analogy to p m , though we 
cannot predict a priori its functional form. Lastly, it is worth 
emphasizing that Eq. ( |B20t is valid to the linear order only 
in metric perturbations as it involves no linearization in the 
matter density fluctuation 8 m , and hence it is valid in principle 
deep in the Newtonian limit, where the validity of the linear 
bias approximation is, however, highly in suspect. 



